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Abstract 

Recent studies show that a negative shock in stock prices will generate more volatil- 
ity than a positive shock of similar magnitude. The aim of this paper is to appraise 
the hypothesis under which the conditional mean and the conditional variance of 
stock returns are asymmetric functions of past information. We compare the results 
for the Portuguese Stock Market Index PSI 20 with six other Stock Market Indices, 
namely the SP 500, FTSE 100, DAX 30, CAC 40, ASE 20, and IBEX 35. In order 
to assess asymmetric volatility we use autoregressive conditional heteroskedasticity 
specifications known as TARCH and EGARCH. We also test for asymmetry after 
controlling for the effect of macroeconomic factors on stock market returns using 
TAR and M-TAR specifications within a VAR framework. Our results show that the 
conditional variance is an asymmetric function of past innovations raising propor- 
tionately more during market declines, a phenomenon known as the leverage effect. 
However, when we control for the effect of changes in macroeconomic variables, we 
find no significant evidence of asymmetric behaviour of the stock market returns. 
There are some signs that the Portuguese Stock Market tends to show somewhat 
less market efficiency than other markets since the effect of the shocks appear to 
take a longer time to dissipate. 

Key words: Asymmetric conditional volatility, stock market returns, threshold 
adjustment, vector autoregression 



Introduction 



It is well known in the financial literature that stock market index returns are 
positively autocorrelated, especially when dealing with high frequency data. 
This fact has been often attributed to non-synchronous trading phenomena 
(see e.g. Fischer [T], Scholes and Williams [2], and Lo and Mackinlay [3]), and 
to time-varying short-term expected returns, or risk premia in the sense of 
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Fama and French [I], and Conrad and Kaul [5]. It is also acknowledged in 
the empirical literature that the unconditional distribution of returns appears 
to be excessively leptokurtic for many stock market index prices. Hamao et 
aZ.[5], and Wei et aZ.[7], for example, found the presence of excess kurtosis in 
the normalized residuals in virtually all mature and emerging stock markets. 
Mandelbrot [8 J , and Fama [9] define excess kurtosis when the return series ex- 
hibit the fat-tailed feature and leptokurtosis when there is observable greater 
excess kurtosis relatively to the normal distribution. Leptokurtosis may pose a 
number of problems when trying to analyse empirically the dynamic behaviour 
of stock markets. In order to deal with these problems, several authors have 
attempted to use more general formulations of the distributions used in empiri- 
cal work (see for example Mandelbrot [Sj, Fama [9], and Nelson [10]). However, 
the subject appears to be far from solved in a totally satisfying manner. 

There are several papers that are studying these topics by applying different 
ideas and techniques ([H];[ll];[l2];[n];[l6];[n];[18]; [19]). 

Another issue obviously related to the previous one concerns the fact that 
short-term stock returns usually exhibit volatility clustering. Links between 
information and volatility of financial assets are of course also central to finan- 
cial economics. Ross [21], for example, argued that volatility could be regarded 
as a measure of the information flow. Moreover, as noted above, empirical re- 
searchers also report significant relationships between flows and volatility. 

Volatility clustering and excess kurtosis may arise in the context of asymmet- 
ric behaviour of price index movements or their corresponding returns. Such 
asymmetry may depend on other economic and financial factors or simply 
be the result of the time-path behaviour of the price index series itself. In 
the latter case, the conditional variance may be an asymmetric function of 
past innovations raising proportionately more during market declines. This 
means that changes in stock prices or returns tend to be negatively related to 
changes in volatility (see e.g. Black [22], and Christie [23]), i.e. the so-called 
leverage effect, where stock price declines increase the financial leverage and, 
consequently, the degree of risk (volatility). 

Authors such as Engle [24], inter alia, defend the idea that financial market 
volatility can be predictable. In order to forecast the volatility, several meth- 
ods have been proposed in the specialised literature, many of which are based 
on regression-like models. Clearly, the classical regression model is unsuited 
when the residuals exhibit features that depart from the basic OLS assump- 
tions. Non-normality, autocorrelation, and heteroskedasticity are some of the 
problems that are typically present in this type of data. In order to deal with 
the problems of autocorrelation and heteroskedasticity, ARCH/GARCH spec- 
ifications are usually well suited. However, they assume that the conditional 
variance has a symmetric behaviour and thus may not fully capture the issues 
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of non-normality. In our context, we shall use TARCH and EG ARCH model 
specifications to assess the extent of non-linear dynamic behaviour underly- 
ing the stock prices data of various markets. The models used in this paper 
belong to special classes of non-linear models that usually generalise the more 
traditional autoregressive conditional heteroskedasticity models. 

Autoregressive conditional heteroskedasticity (ARCH) models have been used 
quite frequently in applied financial research work. Among the most popular 
extensions are the Threshold ARCH (TARCH) model proposed by Zakoian 
[25] . and Glosten et al. [26], and the Exponential GARCH (EGARCH) model 
proposed by Nelson [10J. Engle and Ng [27J indicate that these latter volatility 
models usually outperform standard Generalized ARCH (GARCH) models. 

An alternative way for modelling the volatility of stock market returns is based 
on a Vector Autoregressive (VAR) specification using other variables that may 
help to explain the behaviour of financial markets. In this context we could 
observe e.g. (Sims [|2g|, [2H] and [3U]]). 

A VAR system is a reduced form model that accounts for the effect of lagged 
values of the endogenous variables on the dynamic behaviour of the whole 
system. Each variable is specified as a function of the past values of the other 
(endogenous) variables but does not deppend on their contemporary values. 
VAR models are well suited to deal with situations where an error correction 
mechanism (ECM) can describe the behaviour of the system. In this case a 
corresponding Vector Error Correction model (VECM) can be specified and 
fully estimated using ML procedures. This is important when one needs to 
separate the long-run and short-run components of the dynamic mechanism. 

However, there are several criticisms with respect to VAR models. One is the 
likely loss in efficiency due to lack of apriori information. Another problem 
relates with the large number of parameters to be estimated as the dynamic 
structure of the system increases. This is especially relevant for small data 
sets. Other problems may arise when dealing with impulse response analysys. 

Since stock market returns and other macroeconomic variables may all be 
endogenous in the system, the VAR specification is a quite satisfactory way 
for modelling the volatility. Testing for asymmetries in this context entails 
the use of TAR and M-TAR tests of the VAR residuals (see e.g. Enders and 
Granger [31], and Enders and Siklos [32]). It is important, therefore, to obtain 
comparative results of the relative performance of these two alternative types 
of volatility models, and this is the main aim of the present paper. 
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1 Asymmetry and Volatility Models 

The trade-off between risk and expected return is a foundation of modern 
financial theory, as can be seen in the Capital Asset Pricing Model (CAPM), 
Arbitrage Pricing Model (APT), Portfolio Theory, and Option Pricing The- 
ory. In all these models, volatility plays a central role in determining the fair 
value of a security. To quantify the risk level, several researchers have used 
quantitative statistical indicators, where the standard deviation is one of the 
most popular measures. This statistical measure, however, assumes a symmet- 
ric variation of price levels. Yet, asymmetries are easily observed in time series 
of financial data, and already constitute an important phenomenon. 

For example, there are good reasons to believe that speculative financial asset 
changes have in general an asymmetric behaviour. As mentioned above, the 
leverage effect was found in many empirical studies that analyse the behaviour 
of stock returns. This circumstance highlights the need for using asymmetric 
models when one is analysing data on stock market behaviour. Asymmetric 
behaviour in financial data can be detected with relatively ease, since volatil- 
ity raises more for negative shocks than for positive shock with the same 
amplitude. In order to account for this phenomenon, two extensions of the 
basic GARCH model can be used, among other possibilities: the Threshold 
Autoregressive Conditional Heteroskedasticity (TARCH) and the Exponential 
Generalised Autoregressive Conditional Heteroskedasticity (EGARCH) mod- 
els (see Nelson [10], Zakoian [25], and Glosten et al. [26]). 

Additionally, this study also investigates the relationship between the volatility 
terms for index returns and two macroeconomic variables: Dividend Yield (Dy) 
and the Price Earnings per Share (PER). These variables will be used for TAR 
and M-TAR tests based on a VAR model specification as explained above. This 
approach is similar to that used by Gjerde and Saettem [33] and Rapach [31]. 
Note that the daily frequency of the stock market returns variable imposes 
restrictions on the choice and availability of other macroeconomic variables 
that could also be of relevance but are not available with such a frequency. 



1.1 The Data 

The data used in this study consists of daily stock index closing prices for seven 
stock market indices over a period from January 1993 until December 2002 
(Datastream database). Stock returns were computed as the logarithm of the 
ratio of the index values at time t and time t — 1. We consider the following 
indices: the Standard and Poors 500 (S&P500) to represent the whole US 
economy, the CAC 40 index for France, the FTSE 100 for the UK, the DAX 
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for Germany, the IBEX 35 for Spain, the ASE 20 for Greece and finally, the 
PSI 20 for Portugal. 

The exchanges of Paris, London, Frankfurt, and New York represent four of 
the world's major centres for the trading and distribution of both domestic and 
international equities with special interest for Portuguese investors, and this 
justifies their inclusion in the analysis. Furthermore, France and Germany are 
viewed as "core" European economies, being members of the Exchange Rate 
Mechanism since its inception in 1979, and of the European common currency 
- the Euro - since 1999. The remaining stock markets are more similar to the 
Portuguese stock market and we may expect somewhat similar results for the 
three southern European countries. 



1.2 Methodological Issues and Empirical Results 

The existence of important inter-relationships between the major financial 
stock markets gives support to the use of VAR systems in the context of 
the analysis of volatility. There is the need to understand how shocks and 
volatility in one market are transmitted to other markets. To this end, we need 
to look at the extent to which multi-lateral interaction exists between these 
markets. Thus, after identifying the channels of interaction we can observe the 
implications of volatility co-movements between different markets. 

VAR models have been increasingly used in the analysis of the relationships 
between financial markets. VAR models may be used with both stationary and 
non- stationary data. In the latter case, they provide a convenient framework to 
test for cointegration using the Johansen methodology, but the methodology is 
far more general and is also quite useful when the variables are stationary. On 
the other hand, ARCH/GARCH type models have been widely used to test 
hypotheses concerning the conditional volatility of stock market returns. Both 
methods can be employed to model volatility in the stock market. However, it 
is important to note that employing these different models, each one measures 
and captures different types of volatility. 

TARCH and EGARCH models are especially adequate to model the volatility 
measured in terms of error's variance in the context of asymmetry. On the 
other hand, TAR and M-TAR models systematise the volatility of the depen- 
dent variables in a VAR system, captured from the corresponding residuals. 
The TAR and M-TAR models based on VAR systems are especially designed to 
capture asymmetric multivariate effects in a multi-equation framework, which 
certainly constitutes an important advantage over the single-equation multi- 
variate TARCH and EGARCH models. However, one important disadvantage 
of the former models is that they do not account for conditional heteroskedas- 
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ticity that may be present in the data. 

The first step of our analysis consists of estimating univariate TARCH(1,1) 
and EGARCH(1,1) models in order to assess whether our series exhibit or not 
some type of conditional volatility asymmetric behaviour. For the TARCH 
model we use the following specification of the conditional variance: 

a\ = lo + ae t 2 _i + l4-x<k-i + P°t-i (!) 

where the variance a\ is a function of the past squared residuals e|_!, and of 
its own lagged values o\_ x . The variable dt-i is a dummy variable equal to one 
if e t _i > 0, and equal to zero otherwise, u, a, 7, and f3 are the parameters of 
the conditional variance equation that will be estimated. The equation of the 
conditional variance for the EGARCH model takes the following form: 



log a 2 t =u + 13 log a\_ x + a 



where the symbols described earlier. The results of the univariate 

TARCH and EGARCH models are reported in Table 1. 

As can be seen, the estimates of a are almost all significantly positive, except 
for the US and UK, where we found no significant a coefficients when the 
TARCH model was used. In this model, we obtain good news (e t < 0) and 
bad news {et > 0). This means that the model has differential effects on the 
conditional variance - good news has an impact on a, while bad news has 
an impact on a + 7. We say that the leverage effect exists if 7 7^ and 
an asymmetric effect is observed if 7 > 0. The leverage effect expressed in 
s\_ x d t -\ is associated with the parameter 7 in equation (pQ), and at the same 
time it must reveal statistical significance for both models in Table 1. In the 
EGARCH specification, the asymmetric behaviour exists if 7 < 0. 

Another important result to note is that the estimate of (3 in the EGARCH 
model is less than one for all the markets, which implies that all the moments 
of the relevant statistical distribution exist and that the quasi-maximum like- 
lihood estimators are likely to be consistent and asymptotically normal. Thus, 
relatively to the results reported in Table 1, we may conclude that the leverage 
effect hypothesis cannot be rejected in any case. Likewise, there is evidence of 
an asymmetric behaviour of the stock markets for all the countries analysed. 



e*-i 



+ 7 



£t-i 
<?t-i 



(2) 
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Table 1. TARCH(1,1) and EGARCH(1,1) 



A more detailed analysis of the results reported in Table 1 reveals that the coef- 
ficient a for Portugal is higher than the corresponding values for the remaining 
markets, both for the TARCH and EGARCH models^ These figures imply 
that shocks in the Portuguese stock market tend to have longer durations, or 
higher persistence. As Bala and Prematne [35] state, this circumstance implies 
that the effect of shocks in earlier periods for Portugal tend to linger around 
for a longer period than it does in other stock markets. One possible expla- 
nation for this result is that it may imply that the Portuguese stock market 
shows less market efficiency than the other markets analysed, since the effects 
of shocks take longer time to dissipate. 

1 Except for Greece in the EGARCH model, whose value of a is 0.2602, whereas 
for Portugal is 0.2349. 
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However, one should note that this does not mean that the Portuguese stock 
market is not efficient. Our results should just be analysed in relative terms 
and comparatively with other markets. 

We turn now to consider the results of multivariate TARCH and EGARCH 
models where the Portuguese stock market index is the dependent variable, 
in order to obtain the OLS residuals, and take the remaining market indices 
successively as the independent market index for all possible relations estab- 
lished with the Portuguese stock market. These results may bring new infor- 
mation relatively to the results of Table 1, since it is quite likely that the 
stock markets establish interactions between them. We restrict our analysis 
to the interactions between the Portuguese stock market and the remaining 
markets, successively, and we estimate multivariate TARCH and EGARCH 
models. The results reported in Table. 2 are very similar to those reported 
earlier and lead in general to the same conclusions. The estimated coefficients, 
however, appear now significant in all cases. 

We are now going to test the possibility that an adjustment process occurs for 
different stock market price index combinations with the Portuguese index, 
which have asymmetric cycles in the respective time series. For this test we 
employ the methodology of TAR and M-TAR models, where the indices appear 
as endogenous variables and the macroeconomic factors as control exogenous 
variables. 

The basic starting point for this analysis is the long-run relationship between 
the relevant variables to the model. The long-run equilibrium relationship 
between two time series z t and x t can be estimated as a standard regression 
model z t = a+f3x t +n t , where a and (3 are the estimated parameters, and fi t is 
a disturbance term that may be serially correlated. The parameter f3 gives the 
magnitude of adjustment of z to variations in x, and is the long-run elasticity 
of the two variables if they are measured in logs. If (3 < 1, then shifts in x are 
not fully passed onto z. 

The second step of the methodology focuses on the OLS estimates of p l and 
p 2 in the following error correction model: 



where e t is a white noise disturbance and the residuals from the long-run 
equation are used to estimate Ap t . I t is the Heaviside indicator function such 
that 



A/i t = hp^t-i + (1 - h)p 2 ^t-i + £t 



(3) 



1 if £ t -i > r 
if < r 



(4) 
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Table 2. TARCH and EGARCH multivariate results 



If in (4) £*_ -i — A*t_i) then the model specification illustrated in ([3]) is called 
the threshold autoregressive (TAR) model. It allows for different coefficients of 
positive and negative variations. A sufficient condition for the stationarity of 
p, t is —2 < {pi,p 2 )< 0. This means that the long-run equation is an attractor 
such that n t can be written as an error correction model similar to that given 
in If p x = p 2 then the adjustment is symmetric, which is a special case of 
dSJ) and (jl)). Expression (j3J) can also contain lagged values of Ap. When p t is 
above its long-run equilibrium value, the adjustment is p 1 /i t _ 1 , and if p t _ x is 
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below its long-run equilibrium value, the adjustment is p 2 p t _ v 

If in (jlj) £ t _i = Ap t _ 1 , then the model ([3]) is called the momentum threshold 
autoregressive (M-TAR) model. The M-TAR model allows the decay to depend 
on the previous period change in p, t _ 1 . The value of the threshold r, in our 
case, will be assumed to be zero in all models. 



The TAR model is designed to capture asymmetrically deep movements in 
the series of the deviations from the long-run equilibrium, while the M-TAR 
model is useful to capture the possibility of asymmetrically steep movements 
in the series (Enders and Granger [SI])- For example, in the TAR model 
if —1 < p x < p 2 < 0, then the negative phase of fi t will tend to be more 
persistent than the positive phase. On the other hand, for the M-TAR model, 
if for example \p x \ < \p 2 \ the model exhibits little decay for positive A/i t _ x 
but substantial decay for negative Ap t _ x . This means that increases tend to 
persist but decreases tend to revert quickly toward the attractor. 



Finally, we can perform a number of statistical tests on the estimated coef- 
ficients (and also on the residuals) in order to ascertain the validity of the 
error correction model outlined in (j3J), and subsequently if the adjustment is 
symmetric or not. The relevant tests on the coefficients are H : p x = and 
Ho '■ p 2 = 0, for which we obtain the sample values of the t-statistics; and 
H : p x = p 2 = , for which we obtain the sample values of the F-statistic. 
The restriction that adjustment is symmetric, i.e. Hq : p x = p 2 , can also be 
tested using the usual F-statistic. 



If the variables in the long-run equation are stationary, the usual critical val- 
ues of the t and F distributions can be used to assess the significance level of 
the underlying tests. However, if these variables are integrated of first order, 
one can use the critical values reported by Enders and Siklos [32] to deter- 
mine whether the null hypothesis of no cointegration can be rejected. If the 
alternative hypothesis is accepted, it is possible to test for asymmetric adjust- 
ment using the standard critical values of the F distribution, since p x and p 2 
converge to a multivariate normal distribution (Enders and Granger |31j). 

Table 3 contains the results of the TAR and M-TAR models. As may be seen, 
for the first two tests the null hypothesis is rejected in all cases. Thus, we may 
proceed our analysis with the tests of symmetry and the results are that the 
null hypothesis is not rejected in any of the models considered. 
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Table 3. TAR and M-TAR results 



2 Conclusions 



This paper presents a comparative analysis between four volatility models 
(TARCH, EGARCH, TAR and M-TAR) in order to ascertain their power 
to capture asymmetric cycles in the underlying time series. Our results indi- 
cate that while the TAR and M-TAR models do not identify the presence of 
asymmetry in any market, the TARCH and EGARCH models have captured 
asymmetry in all cases under analysis. 
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Another relevant conclusion from the multivariate TARCH and EGARCH 
results is that the Portuguese stock market appears to reveal higher a values 
than the remaining markets, for both models. These figures imply that shocks 
in the Portuguese market, in the later periods, tend to have longer duration 
periods. This circumstance implies that the effects of shocks in earlier periods 
tend to linger around for a longer period than it does in other stock markets. 
This means that for the Portuguese stock market the effects of shocks take a 
longer time to dissipate. 
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